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Abstract 

Vibrational energies and wave functions of the triplet state of the ion have been determined. 
In the calculations, the ground and first excited triplet electronic states are included as well as the 
non-Born-Oppenheimer coupling between them. A diabatization procedure transforming the two 
adiabatic ab initio potential energy surfaces of the triplet-H^ state into a 2 x 2 matrix is employed. 
The diabatization takes into account the non-Born-Oppenheimer coupling and the effect of the 
geometrical phase due to the conical intersection between the two adiabatic potential surfaces. 
The results are compared to the calculation involving only the lowest adiabatic potential energy 
surface of the triplet-H^ ion and neglecting the geometrical phase. The energy difference between 
results with and without the non-adiabatic coupling and the geometrical phase is about a wave 
number for the lowest vibrational levels. 
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I. INTRODUCTION 


Hydrogen is the most abundant element in the universe. The most important hydrogenic 
compounds are the dihyrogen molecule and its protonated form, H^“. The latter is a strong 
acid and very reactive, but does exist in large amounts in low density environments such 


as the in 
Oka 


1C 111 

q a 


ierstellar space. The role it plays in interstellar chemistry has been reviewed by 


2]. The ground electronic state of is a singlet, X 3 Ai, with the minimum of the 
potential energy surface (PES) at the equilateral geometry (EG) of three nuclei with the 
internuclear distance r = 1.6504 bolir (do below). 

Above the ground singlet state of H 3 , there is a metastable triplet state, d 3 E+. It has 
never been observed but has been predicted to exist in 1974 by Schaad and Hicks [sj] and 
explored theoretically. The Erst complete potential energy surfaces (PES) and vibrational 
and rovibrational calculations were published only 27 years later and independently by 
two groups 4, 5]. The lowest triplet PES is characterized by three minima at symmetric 
linear configurations with rq = r 2 = r 3 /2 = 2.4537 do- They are separated by barriers of 
2596 cm” 1 , with configurations rq = 1.99 a 0 ,r 2 = r 3 = 5.336 a 0 corresponding to T-shaped 
complexes of a H^“ ion and a distant hydrogen atom. The three dissociation channels to 
H+ + H are at 2946.8 cm 1 . At EG, the state is degenerate with the second triplet state 
for symmetry reasons. The curve of degeneracy has a minimum at E = 14934 cm ” 1 above 
the dissociation energy or E = 17880.8 cm ” 1 above the minimum, with internuclar distance 
rq = rq = r 3 = 3.610 do- As the molecule departs from the equilateral configuration, 
the degeneracy is lifted according to the Jahn-Teller theorem jl5]. Therefore, the two triplet 
electronic states could be viewed as two components of the 3 E' electronic state in the D 3h (M ) 
molecular symmetry group, which intersect conically at equilateral configurations. For the 
simplicity of discussion, we will be referring the two lowest electronic triplet states as the 


lower (/) and upper ( u ) states. Individual fits o: 
rovibrational calculations have been published 


’ either of the two states and corresponding 
7-9]. A two-surface potential has also been 


reported 10]. For an overview on triplet H 3 see Alijah and Varandas fll ]. 

Vibrational energies obtained taking into account only the lowest triplet state (one-state 
calculations) and neglecting the interaction with the upper state may not be precise even 
though the upper state is located at high energy, in the continuum of the lower state. In this 
study we discuss the effect of the coupling with the upper state on the vibrational energies 
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of the lower state and in particular the effect of the geometrical phase. 


II. THEORETICAL APPROACH 


A. Diabatization representation of the potential 


The triplet-Hjj" PES calculated by Varandas et al. [uj is utilized in the present work. It 
provides an accurate representation of the two components, in which their degeneracy along 
the line of the conical intersection and the linear splitting in energy with the symmetry 
breaking vibrational coordinate near the intersection seam is guaranteed by construction. 
This PES is shown in Fig. [1] In order to take into account the non-adiabatic interaction 
between the two components, we employ an approach similar to the one discussed in Ref. 12]. 
Namely, we transform the two-component adiabatic (Born-Oppenheimer approximation) 


PES of Ref 


'• M 


into a diabatic form, where the potential is written in the basis of two 


degenerate electronic functions |+) and |—). For the diabatic electronic basis, we use the 


electronic states evaluated at the EG 


13]. At the EG, the electronic functions of the triplet- 


ion are doubly degenerate so we can choose any linear combinations of the two functions 
as the basis. Following Longuet-Higgins [13J], we choose the linear combinations that are 
orthogonal to each other and transform in the C 3v (M) symmetry group as E + and E_ 
components of the E representation: 


(123)£+ = uE + ; (12 )E+ = E _ ; 

(123)E_ = u*E_ ; (12 )E_ = E + , (1) 

where (123) and (12) are the cyclic and binary permutations of H^“ nuclei [m)] and u = e 27r */' 3 . 
In this chapter, we use the symmetry labels of the molecular symmetry group C 3v (M), since 
the odd-parity representations of D 3h (M ) are not needed. 

In the diabatic basis, the potential of interaction between the three nuclei is written 
as a 2 x 2 matrix with matrix elements dependent on geometry. Because the energy of 
two components of the potential is the same at EG, the diagonal elements of the diabatic 
potential are the same. In a general form the potential is written as 


( A Ce^ 
y Ce-V A 


( 2 ) 


3 





FIG. 1: Born-Oppenheimer (adiabatic) potential energy surfaces of the two components of the H 3 
triplet state shown as a function of two asymmetric stretch coordinates, Q x and Q y of Eq. (J3]) for 
a fixed value of the symmetric stretch coordinate Q± = 5.3 ao- The three minima of the shown 
part of the lowest PES are situated at linear isosceles geometries of the three nuclei. 
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where A, C, and / are real-valued functions of internal coordinates. The functions A and C 
are transformed in the C 3v (M) molecular symmetry group according to the A\ irreducible 
representation. The phase / of the coupling has the following property under the C 3 sym¬ 
metry operator of the C 3v (M ) group: C 3 f = f + 2n/3 , i.e. the non-diagonal elements 
transform as E + and E_. The A and C functions are uniquely determined from the two 
Born-Oppenheimer components, Vi and V u , of the triplet state PES: A = (y + 14)/2 and 
C = (Vi — V u )/2. The actual form of the function / can be derived near the EG using the 
diabatization procedure suggested by Longuet-Higgins [lsj]. I 11 the diabatization procedure, 
the phase of the non-diagonal coupling is equal to the phase a of asymmetric distortion 
jl3 . 15] in symmetry-adapted coordinates. The symmetry-adapted coordinates Qi,Q x and 
Q y are related to three internuclear distances r 1 , r 2 ,r 3 as 


Qi = (ri + r 2 + r 3 )/V 3, 

Qx = (2r 3 - r v ^ r 2 )/v / 3, (3) 

Q y = r 2 -r 1 . 


Q x and Q y become degenerate at D 3h configurations and may be expressed more conveniently 
in polar form, 

P = \J Ql. + Q't 5 

Q x — P cos a , (4) 

Q y = p sin a . 


As suggested by Longuet-Higgins, we set f = a. Although this form of / is derived near the 
EG, we will also use it in the region far from the EG. It is a good approximation because 
( 1 ) the transition between adiabatic states occur only near the EG, where energies of the 
two PES components approaches, and (2) far from the EG, the phase factor e ±l ^ does not 
play a significant role as long as the symmetry of / in the group C 3v (M), mentioned above, 
is satisfied. 

Notice that the basis defined above is diabatic with respect to the asymmetric stretch 
coordinates only but it still depends on the symmetric stretch, i.e it is adiabatic along the 
coordinate. However, for the simplicity of discussion, we refer this basis as diabatic, even if 
it is only partially diabatic. 
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The solution of the Schrodinger equation for the vibrational motion with the potential 
of Eq. (J2J) is a function T having two components - 0 ± corresponding to the two diabatic 
electronic states. In the adiabatic basis, corresponding to the |Z) and | u) electronic states, 
the two components ipi >u of X I> have the form 

0; = (0_e* //2 + 0 + e _i//2 )/\/2, 

ip u = i(0_e*^ 2 — 0+e - *-^ 2 )/ V2 . (5) 

The C 3 symmetry operator applied three times to the molecule returns the three nuclei 
back to their original relative configuration. However the components change sign 
—> —^i jU under the Cf operator because /—>■/ + 2tt. It is a well-known property of 
adiabatic states when there is a conical intersection. The property is often referred to as 
geometrical or Berry phase effect. Because the adiabatic electronic wave functions |Z) and 
| u) also change sign, the total vibronic wave function 

'k = ipi\l) +i> u \u) = -0+1+) + V’-l-) (6) 

is unchanged after the Cf operation. 

Therefore, the components are two-valued. This complicates the solution of the 
vibrational Schrodinger equation in the adiabatic basis. I 11 this situation, the diabatic basis 
is preferable because the diabatic components are single-valued. In addition, in the diabatic 
representation, there is no need to deal with first and second derivatives of electronic wave 
functions in order to account for non-adiabatic couplings between the two adiabatic PES. 


B. Nuclear dynamics 


The hyperspherical coordinates R , 9 ,0 are specially adapted to solve the vibrational 
Schrodinger equation for Hjj“. Solving the Schrodinger equation, we use the two-step ap¬ 
proach 16]. First, we solve the equation for a fixed hyper-radius R, 3 


H Rj (p j: a(9, 0) = UaiRj^aie, 0) , (7) 

where H R . is the two-channel vibrational Hamiltonian with the hyper-radius fixed at Rj. The 
solution of the above eigenvalue problem gives eigenstates <fj t a(9, 0 ) depending on the two 
hyperangles, which are usually called hyperspherical adiabatic (HSA) states, and eigenvalues 
U a (Rj), called HSA potentials. 
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Once the HSA states are determined, the complete vibrational wave function 9, 0) 
of the two-component triplet potential is represented as an expansion -0(77, 0,0) = 
J2j, a yjA R i 0) c ja in tlie basis of following functions y ha = 'K j (R)cp jta (9, 0), where <p jia (9, 0) 
are calculated at the hrst step, tt j(R) are DVR basis functions along the hyper-radius, local¬ 
ized at values Rj of the hyper-radius. Coefficients Cj a are unknown. The expansion functions 
yj t a are not necessarily orthogonal to each other because HSA states tpj, a (9, 0) evaluated at 
different Rj are not orthogonal. Therefore, in the basis y ]a , the vibrational Schrodinger 
equation for 0(i7, 9, 0) takes the form of a generalized eigenvalue problem for coefficients Cj a 


F 1 

h 2 d 2 

[(^\ 

2/i dR? 


Cj'a' 


) Oja,j'a' “t - I ^a(-R) 

= E n ^ (tTj I TTj') OjaJ'a 1 Cj'g' , 


( 8 ) 


where fi is the three-body reduced mass of the ion; Oj a j ' a 1 are overlaps between HSA states 


C> 






(9) 


III. BOUND STATES OF THE H+ TRIPLET STATE 

At equilibrium, the H0 molecule in the triplet state is linear. Three equivalent min¬ 
ima exist due to permutationial symmetry of identical nuclei, separated by barriers. Since 
the barriers are not too high, tunneling is a feasible operation. Furthermore, most of the 
vibrational states are located above the barriers. Symmetry classification of the rotation- 
vibrational states of this floppy molecule then requires the use of the full three-particle 
permutation-inversion group S 3 x /, which is isomorphic with the /^ (M) molecular sym¬ 
metry group. Pure vibrational states ( J = 0) must have even parity and hence transform as 
one of the three irreducible representations of the subgroup C 3v (M), A\ , A! 2 , and E'. Below, 
we will omit the primes. 

To obtain vibrational energies and wave functions, hrst, we have obtained the HSA wave 
functions and energies. In the numerical calculations, we used B-splines [itJ to solve the 
hyperangular part of the Schrddinger equation. I 11 a typical calculation, we used 180 B-splinc 
intervals along angle 9 , which changes in the interval [— 7T, 7t] . The angle 0 changes in the 
interval [0, 27t] . However, because of the symmetry of the potential as a function of 0, each 
HSA state ip (as well as the total vibrational functions T) belongs to one of the three irreps 
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Ai, A 2 , and E and, as a result, has a certain periodicity as a function of 0. The periodicity 
allows us to reduce the [0, 27r] interval: To obtain HSA states of the A\ and A 2 irreps, the 
interval [—7 t/2 ,7 t/6 ] is sufficient, and for states corresponding to the E irrep, the interval 
[—7 t/2 ,7 t/2 ] is enough, provided that appropriate boundary conditions are applied to the 
HSA states at ends of the intervals. For the Ai states and one ( E a ) of the two components 
of the E irrep, derivative of wave functions with respect to 0 are zero at the ends. For the 
A 2 states and the second component (F0) of the E irrep, the wave functions itself are zero 
at the ends. In a typical calculation, we use 100 B-splines along the 0 to obtain A\ and A 2 
states and 211 B-splines to obtain E states. 

The wave functions and have two electronic components. By def¬ 

inition, the electronic diabatic states |+) and |—) transform as components E + and E_ 
of the E irrep. Therefore, the irrep T v n, r of the vibronic function is related to the ir¬ 
rep T vi b of the vibrational components i\)± as T vibr = ^mb x E. The vibrational adiabatic 
components < Pj >a ,e{9 ,0) with e = l or u exhibit the property of geometrical phase: If the 
geometry of three nuclei is changing continuously in such a way that it makes a closed loop 
in the space of geometries and the equilateral configuration stays inside the loop, the wave 
function changes sign when the geometry returns to the initial one. Therefore, adiabatic 
components of total wave functions could be considered as two-valued: An arbitrary config¬ 
uration of the Hjj“ nuclei described by an angle a of distortion (or by a hyperangle 0 ) can 
equally be described by a + 27r (or 0 + 2n). Each adiabatic component is then two-valued 
<Pj, a ,e{8i 0) = — ( fj,a,e{9,(p + 27T). Because the electronic adiabatic states are also two-valued, 
the total wave function stays single-valued. 

Examples of diabatic and adiabatic vibrational components of <pj ta (9,(j)) are shown in 
Figs. [2] and [3] for a fixed value of R — 4.7 ao for T vi b = A\. Figures [4] and [5] shows several 
HSA states for T vib = E. Wave functions of bound states in the diabatic representation of 
complex-values electronic states |+) and |—) are in general also complex. However, adiabatic 
states |Z) and | u) are read and, therefore, the overall phase of the adiabatic vibrational 
components could be made real. In Figs. [2] and 0] we show real and imaginary parts of the 
two components of lowest three (a = 1,2,3) HSA (pj >a (9,<f) in the diabatic representation. 
Figures [3] and [5] show only real part of the two components of lowest six (a = 1, • • • , 6 ) HSA. 
The imaginary part of the adiabatic components is zero. 

The property geometrical phase is demonstrated by vibrational adiabatic components of 
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FIG. 2: The first three (a = 1,2,3) HSA wave functions ^>(0, </>) of the A\ irreducible representa¬ 
tion of H^" in the diabatic representation obtained for a fixed value of hyper-radius Rj = 4.7ao- Each 
HSA wave function has two complex-valued components, 99 + and 99 _ similarly to two functions ?/>+ 
and 9 !of Eq. ([ 6 ]). They are shown as a column of four graphs for each a. 

r vib = Ai and E in Figs. [3] and [5j Only one sheet of the two-valued functions is shown. 
The equilateral geometry is the branch point for the cut. The configuration space of two 
hyperangle is cut from the branch point straight down to the linear geometry along the line 
<f> = —7 t/2 or a = tt. The choice of the cut is, of cause, arbitrary as long as it goes from the 
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FIG. 3: The first six (a = 1, • • • ,6) Ai HSA wave functions </>) shown here in the adiabatic 

representation. The first three functions are the same as in Fig. El In the adiabatic representation, 
the wave functions could be made real and, therefore, the imaginary parts are zero and are not 
shown here. Notice that the two adiabatic components of the vibrational wave functions change sign 
when the hyperspherical angle (j) goes through value —7r/2, which is in agreement with to Eq. (f5|). 
Also, notice that the white color corresponds to (fj^ a = 0 in all graphs, but the correspondence of 
intensity of red and blue colors to actual values of (pj, a is slightly different in the graphs. 
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FIG. 4: The first three (a = 1,2,3) HSA wave functions (pj >a (6,(j)) of the E a irreducible represen¬ 
tation in the diabatic representation calculated for Rj = 4.7ao- Each HSA wave function has two 
complex-valued components, 99 + and similarly to two functions ?/>+ and ip— of Eq. (0. They 
are shown as a column of four graphs for each a. 

branch point to a linear geometry. As it is demonstrated by Figs. [3] and El the adiabatic 
components indeed change sign at the cut. In Fig. [5] it is less evident, except a = 6 , because 
the wave functions with a = 1, • • • , 5 are almost zero near the cut. 

To assess the effect of the coupling with the upper PES, we have made two types of 
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FIG. 5: The first six (a = 1, ■ ■ ■ , 6 ) E a HSA wave functions Lpj^ a {0, <p) shown here in the adiabatic 
representation. The first three functions are the same as in Fig. [2 

calculations. First, we determined HSA states and vibrational energies using only the lowest 
adiabatic potential Vi - the one-state calculation. In the second type of calculation, we solve 
the Schrodinger equation for the coupled potential of Eq. ([2]) - the two-state calculation. The 
difference between energies of obtained vibrational states reflects the effect of the coupling 
between the two adiabatic surfaces on the vibrational energies of Vj. Figure [ 6 ] shows the HSA 
energies U a (R) obtained in one- and two-state calculations. The HSA potentials curves are 
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Hyper-radius ( a Q ) 

FIG. 6: Hyperspherical adiabatic potentials of A\ vibrational symmetry obtained for the triplet 
state of Hjjb Due to the fact that the conical intersection between the surfaces of the two compo¬ 
nents of the triplet state is situated at a relatively high energy compared to the minimum of the 
potential, the difference between the HSA potentials obtained in one- and two-state calculations is 
small, about 3 x 10~ 5 E\ u and cannot be seen at the figure. 

very similar to each other in one- and two-state calculations in the energy region, shown in 
the figure. For example, near the minimum of the lowest HSA curve, the energy difference 
between the two calculations is about 3 x ICE 5 E h . 

Once the HSA energies and wave functions are calculated, we solve the hyper-radial 
part (Eq. [9]) of the Schrodinger equation. Energies of obtained vibrational levels in one- 
and two-state calculations are listed in Table HI In the table, symbols of irreducible rep¬ 
resentations refer to symmetry of the vibrational part in one-state calculation, and to the 
symmetry of the diabatic vibrational components in the two-state calculation. In addition, 
we provide an identification of the vibrational states in terms of D^^M) linear molecule 
spectroscopic quantum numbers (iq, v^, u 3 ), where V\ and v 3 correspond to the symmetric 
and antisymmetric stretching vibrations, and v 2 to the degenerate bending vibration, with 
l the associated vibrational angular momentum. Such a classification is not exact, since 
the molecule is not localized in a single potential well. Rather, linear combinations of the 
three equivalent vibrational wavefunctions must be formed, which then give rise to a one- 
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dimensional and a two-dimensional representation in the C$ V (M ) group. For not too highly 
excited vibrational states, the linear molecule quantum numbers could be useful to label the 
vibrational states if supplemented by a degeneracy index (A\/A 2 or E). A discussion of the 


Results of the one-state calculation may be 
ldl ]. where the same surface Vj was used. The 
could not be converged to 


quantum numbers can be found in Refs. % 
compared with that a previous calculation 
results agree well except for the highest states, which in Ref. jlo| 
the same accuracy as the lower states. 

Comparing now the relative energies of the one-state and two-state calculations, they 
agree to within fa lcm -1 for levels up to fa 2300 cm -1 . These small differences can be 
attributed to the neglect of non-adiabatic effects in the one-state calculation. It is worth 
noting that only the first three states are located below the electronic energy of the H^“ + H 
channel (E = 1228.5 cm -1 ), while the large remainder lies between this threshold and the 
vibrational zero point energy of H^“. 

The splitting in energy between the one-dimensional and the two-dimensional components 
of a linear molecule type vibrational state is expected to be small for the lowest energy levels 
because the tunneling is weak. For more excited vibrational levels the splitting is expected to 
increase, which is indeed observed in one- and two-state calculations as evident from Table U 
The tunneling between the three potential wells is more effective in the two-state calculations. 
It happens because the non-adiabatic coupling with the upper potential surface V u makes 
tunneling easier if compared to the calculation involving only the Vj surface. For example, 
the lowest excited state (0, 0°, 0 )E in one-state calculation is just about 2 x 10 -4 cm -1 above 
the ground vibrational level (0, 0°,0)Ai. In the two-state calculation, the energy difference 
between (0, 0°, 0 )E and (0, 0°, 0)Ax is 1.04 cm -1 . For the second state, (0, 0°, 1), the splitting 
is 0.05 cm -1 in the two-state calculation. The higher vibrational states are already located 
above the barriers. 

Not all of the computed states shown in Table U are allowed for this system with three 
identical protons. Since protons are fermions, the total wavefunction must be antisymmetric 
under permutation of the position of two nuclei and hence transform as A' 2 or A 2 in the 
group i(M). The total wavefunction T tot has contributions from the nuclear spin (ns), 
electronic spin ( es ) and the spatial parts, electronic (el), vibrational (vib) and rotational 
(rot). If rotation is neglected (J = 0), classification in the C^ V (M) symmetry group is 
sufficient. The electronic spin function is totally symmetric and transforms as A 1 . It suffices 
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TABLE I: Energies of vibrational levels obtained in one- and two-state calculations. The origin is chosen at the energy of the ground 
vibrational level (0,0°, 0). In one-state calculation, the zero-point energy is Eq = 1722.043cm -1 , the energy of H^~(u = 0 ,j = 0)+H 
dissociation is E t hres = 2382.0 cm -1 above (0,0°, 0). 
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to analyse the nuclear spin, electronic and vibrational symmetries. 

Two nuclear spin functions may be formed, one having 7 = 3/2 (ortho spin) and the 
other 7 = 1/2 (para spin). They transform as A\ and E, respectively. The symmetry of 
the electronic functions in the adiabatic representation is A 2 for the lower surface and Ai 
for the upper, while in the diabatic representation they transform as two components, E + 
and E_ (or E a and 7/,), of the degenerate representation E. Hence, with respect to their 
symmetry under the (12) permutation operation, the lower adiabatic state is antisymmetric 
and the upper symmetric. This is different for the diabatic representation, where the two 
components are transformed between each other, see Eq. ([TJ) . 

In the adiabatic representation, we have for ortho nuclear spin the direct product r e * x 
T ns = A 2 x Ai = A 2 which requires the vibrational state to have A 1 symmetry so that the 
total function transfroms as 


Ttoi — T e i x T ns x T vi b — A 2 x Ai x Ai — A 2 (10) 

For para spin, T e ; x T ns = A 2 x E = E, which is compatible with E vibrational symmetry, 
since 

Ttot — T e z x T ns x T v ib = A 2 x E x E = A 2 (+Ai + E) , (11) 

where the Ai and E irreps are taken in parentheses because they are not allowed for the 
complete wave functions of Hj/. 

We now turn to the diabatic representation, where we have T e / = E. Hence for ortho 
nuclear spin T e ; x T n , s = E x A\ = E which selects T,,^ = E, since 


f(oi - fei x T ns x T v ib — E x A\ x E — A 2 (+Hi + E) 


( 12 ) 


For para nuclear spin we obtain T e ; x T ns = E x E = Ai + A 2 + E which may be combined 
with any vibrational symmetry, 

( 


v tot = r pi x r M x r „ ih = 


el ^ ns 


vib 


E x E x A 1 = A 2 (+Ai + E) 
E x E x A 2 = A 2 (+A 1 + E) 
E x E x E = A 2 (+Hi + 3 E) 


(13) 


An important consequence of the geometrical phase is therefore its alteration of the spin 
statistical weights of the vibrational functions which we present in Table [TT1 


16 




TABLE II: Nuclear spin statistical weights of the vibrational wave functions obtained in a one-state 
and a two-state calculation. The geometrical phase present in the two-state calculation alters the 
spin statistical weights. In particular, it makes A -2 states allowed. This difference could be crucial 
to diagnose the geometrical phase experimentally. 


F vib 

W(lst.) 

1 W(2st.) 

At 

4 

2 

A 2 

0 

2 

E 

2 

4 


IV. CONCLUSION 

Vibrational energies and wave functions of the triplet state of the ion have been 
determined. In the calculations, the ground and first excited triplet electronic states are 
included as well as the non-adiabatic coupling between them. The calculations have been 
performed in a diabatic basis using a diabatization procedure suggested by Longuet-Higgins 
transforming two adiabatic ab initio potential energy surfaces of the triplet-Hg state into a 
2x2 matrix. The suggested theoretical approach allows us to account for the non-adiabatic 
coupling between the ground and the first excited electronic states of triplet-Hg . It also 
includes the effect of the geometrical phase due to the conical intersection between the 
two adiabatic potential surfaces. The results are compared to the calculation involving 
only the lowest adiabatic potential energy surface of the triplet-H^ ion and neglecting the 
geometrical phase. The energy difference between results with and without the non-adiabatic 
coupling and the geometrical phase is about one wave number for the lowest vibrational 
levels. However, in the presence of the geometrical phase the spin statistical weights of the 
vibrational states are altered. These findings are consistent with the expected small effect 
on the vibrational energy levels as the excited triplet state is located far in the continuum 
of the lower state. They are also consistent with the fact that the geometrical phase is 
topological and therefore manifests itself even if the two electronic states are well separated 
in energy. 
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